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We consider a very simple toy model for a spatially varying 'cosmological constant', where we 
are inside a spherical bubble (with a given set of cosmological parameters) that is surrounded by a 
larger region where these parameters are different. This model includes essential features of more 
' realistic scenarios with a minimum number of parameters. We calculate the luminosity distance 

in the presence of spatial variations of the vacuum energy density using linear perturbation theory 
and discuss the use of type la supernovae to impose constraints on this type of models. We find 
that presently available observations are only constraining at very low redshifts, but also provide 
independent confirmation that the high red-shift supernovae data does prefer a relatively large 
positive cosmological constant. 
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The Standard Cosmological Principle (SCP) states that the universe is spatially homogeneous and isotropic on 
large scales. It is one of the cornerstones of modern cosmology, in the sense that if it were not true, then the standard 
results for the most basic properties of the universe (such as its geometry, content, or age) would not hold and would 
Oh | have to be re-worked from scratch, quite possibly suffering rather drastic changes. There is some supporting evidence 
q ■ for the SCP on scales very close to the horizon though it can hardly be called definitive. 

However, the theoretical situation is not totally unambiguous either. On one hand, it has been claimed Q that the 
SCP doesn't follow from the Cosmic Microwave Background Radiation (CMBR) (near) isotropy and the Copernican 
Principle, as is usually assumed. This is important because in that case it follows that the SCP cannot be assumed 
J> \ to hold based on the presently available observational data. On the other hand, it is known that there are ways |^]] 
in which the SCP could be evaded that would be difficult, or even impossible, for us to notice at the present time. 

A simple example which the present authors have considered in the past is that of a late-time phase transition 
producing wall-like defects which separate regions with different values of the cosmological parameters, notably the 
matter and vacuum energy densities (and hence also the Hubble constant). Since in such models the differences in the 
cosmological parameter values in the different regions tend to increase with cosmic time after the phase transition that 
produced them, and since cosmological observations necessarily look back at earlier times, detecting such differences 
is not as simple as one might expect, and hence constraints on these models turn out to be surprisingly mild. 

In addition, cosmological observations in such models will be faced with the same type of limitations as were pointed 
out in H for the case of quintessence models (see also ppp[| for different perspectives). It should be noticed that 
although the two types of models have rather different motivations, they will be rather similar observationally, since 
in both cases one is effectively dealing with an equation state of the universe which varies as a function of the redshift. 
However, an important additional feature which needs to be taken into account in models with different domains is 
the temperature jump due to the relative velocity of comoving observers on either side of the domain wall. Also, 
these models will not be isotropic and so, in general, results will vary as a function of the direction on the sky. This 
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means that problems with systematic errors due to supernovae evolution or dust are less severe than in the context 
of quintessence models. 

In this paper we will use the recent measurements of the luminosity redshift relation using supernovae out to z ~ 1 
to constrain this type of models. We should emphasise at the outset that we will only be dealing with a very simple 
'toy model', but we nevertheless believe that it still includes crucial features of more realistic scenarios with a minimum 
number of parameters. 

Supernovae have long been recognised as a crucial step in the cosmological distance ladder, and hence as a useful 
tool to estimate cosmological parameters via the Hubble diagram. Recent observations of type la supernovae by two 
independent teams PHl^t (the Supernovae Cosmology Project and the High-z Supernova Search Team) out to redshift 
z ~ 1 provided, together with CMB data, some fairly strong evidence for a recently started phase of acceleration in 
the local universe, with the preferred vacuum and matter densities being J1a ~ 0.7 and f2 m ~ 0.3. This conclusion 
is based on the observed faintness of high-redshift supernovae, relative to their expected brightness in a 'standard' 
decelerating universe. It should be emphasised that these measurements are local and can not be extrapolated all the 
way to the horizon. Therefore they do not, on their own, imply that we have entered an inflationary phase 

In the following section we introduce our toy model and discuss its possible shortcomings. We then study the 
evolution of linear perturbations in our model in Sect. and derive an expression for the luminosity distance as a 
function of the red-shift in Sect. IV. Finally, Sect. ^ contains a thorough discussion of our results, and we conclude 



in Sect. VI 



II. THE MODEL 

In a recent article j?J we have shown that large sub-horizon inhomogeneities may be generated if a network of domain 
walls permeates the universe, dividing it in domains with slightly different values of the vacuum energy density and 
other cosmological parameters. The typical size of these regions is determined by the dynamics of the network of 
domain walls and is expected to be close to the horizon scale. The necessary condition in order for the model to be 
observationally viable is that the domain walls are formed in a late-time phase transition. 

This condition is required for two different (though related) reasons. Firstly, the cosmological parameters will be 
different in the different domains, and the differences tend to increase with time, so they must not be so big as to 
make it observationally obvious at the present time. Secondly, the domain walls can themselves be cosmologically 
important (or even disastrous), and in order to avoid this one requires that their energy scale is sufficiently low so 
that they do not contribute in a significant manner to the CMB anisotropics. 

Here we study a simplified model in which the universe is made up of a spherically symmetric region (domain) , which 
is surrounded by another region with a different vacuum energy density (which we will call V- and V+ respectively). 
We assume that we live in the centre of the inner region. Moreover, one assumes that the thin region separating the 
two domains considered (domain wall) does not generate relevant CMB fluctuations. This happens if the potential of 
the field is small enough at the origin. We also require that the domain walls have no non-trivial dynamics, which is 
a good approximation if friction is important [p~5|— [l7|| . 

The vacuum density will be parametrised by 

Oa^, (1) 

Pc 

where p c is the critical density, and we define AOa as 

A fi A (r) = ^ = P^- V -. (2 ) 

Pc Pc 

where p\{r) is the vacuum energy density at the point in question. Hence, this can have two possible values: if we 
are inside the inner region, and (V+ — V-)/p c in the outer domain. 

We emphasise that this is a highly simplified model. In more realistic models, domain walls will have a non-trivial 
dynamics and the domains will have many different shapes. Furthermore we are also not expected to be exactly at 
the centre of a given domain. In general, such cases would have to be dealt with numerically. We will analyze some 
of these effects elsewhere However, this simplified model still incorporates some of the crucial features of more 
realistic models with a minimum number of new parameters (the red-shift of the domain wall z* and the difference 
between the two vacuum energy densities ASIa). In the next sections we will show how high-redshift supernovae 
can be used to constrain combinations of these parameters. We shall also discuss the validity of these results in the 
context of more realistic models in which a network of domain walls is present. 
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III. EVOLUTION OF COSMOLOGICAL PERTURBATIONS 



In the conformal-Newtonian gauge, the line-element for a flat Friedmann-Robertson- Walker background and scalar 
metric perturbations can be written as 

ds 2 = a 2 {ri) [(1 + 2$)c 2 dr/ 2 - (1 - 2<$>){dr 2 + r 2 d0 2 + r 2 sin 2 Od<t> 2 )} , (3) 

assuming that the anisotropic stresses are small. Here, $ is the metric perturbation, c is the speed of light in vacuum, 
a is the scale factor, r\ is the conformal time, and r, 9 and 4> & r e spatial coordinates. 

Given that the vacuum energy becomes dominant only for recent epochs we shall be concerned with the evolution 
of perturbations only in the matter-dominated era, neglecting the contribution of the radiation component. The 
evolution of the scale factor a is governed by the Friedmann equation^ 

H 2 = fl^a- 1 + fl° A a 2 . (4) 

Note that the background matter and vacuum energy densities at an arbitrary epoch can be written as 

nm = ns, +7$ fl a (5) 

and f^A = 1 — Q m - 

The linear evolution equation of the scalar perturbations is given by (see for example fll9|| ) 

$ + 3H<S> + [2H + H 2 ]$ = 4:TrGa 2 5p = -^a 2 An A , (6) 
where dp is the pressure perturbation and 

2H = -fl^a -1 + 2n° A a 2 . (7) 

Given that the source term in the outer domain (oc a 2 A£ljy) is only important near the present time, we shall assume 
the following initial conditions for eq. (^|): 

$(0) = 0, $(0) = 0, (8) 

both in the inner and the outer regions. We note that given that the source term is absent in the inner region, the 
metric perturbation is always zero there. 

The density perturbation, <5, in the outer domain is given by Jl9| 

j=^ =2 ^ = 2 w -a [v a $ _3 W ^_3 W a $]) (9) 
p c no 

which simplifies to 

5/2=^ = -(H- 1 i> + <S>), (10) 

given that V 2 $ = except at the domain walls. The relationship between S and the fractional perturbation in 
the expansion rate follows directly from the Friedmann equation. It is straightforward to show that if the metric 
perturbations are small the outer domain behaves as having an effective f2 A ff and fi^f given by 

= (l - ("a + A^ A ) (11) 

and similarly fi^f = 1 - ilf. 



A dot denotes a derivative with respect to conformal time, TL = a/a, the index '0' means that the quantities are to be 
evaluated at the present time and we have taken ao = 1 and TCo — 1 (so that the conformal time is measured in units of H-q 1 )- 
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IV. THE LUMINOSITY DISTANCE 



We can now proceed by evaluating the luminosity distance relation in the presence of spatial variations of the 
cosmological parameters. Recall that the luminosity distance to a given source is given by: 

4 = -Af (12) 

where C is the luminosity of the source and T the measured flux. It follows directly from the perturbed flat FRW 
metric given in eqn. (|^) that, to first order in the metric perturbation, the luminosity distance, g?l, to an object with 
comoving coordinate r — r\ at a red-shift z is given by: 



d L (z) ~ (i + Z )[i - HmMn ~ (i + z)[i - $(%,o)]c / [i + 2§tf, \m - v»W, (13) 

J V (z) 



where <&(r)',r') = $(77 = rf , r = r'), 770 is the conformal time at the present time, n defines an arbitrary direction on 
the sky and the observer is assumed to be at r = 0. 

Analogous corrections will also arise for the Sachs- Wolfe effect. Here, in the presence of the scalar metric pertur- 
bations defined by eqn. m) there is an additional shift in the temperature of the source given by: 



AT _ f'° 



rVo 

/ dry$(?7, |?7o-r/|n) + 2$(?7o,0)-2$(?7 e ,|77o-77 e |n). (14) 
J n« 



where r\ e is the conformal time when the light was emitted. Within a given domain AT/T is obviously zero. However, 
if the source is in the outer domain there will be a temperature jump at the domain wall with 

AT 

— = -2$+. (15) 

where <f> + is the value of <f> in the outer region at the time when the light crossed the domain wall (here we are 
assuming that $ = in the inner region). This temperature shift is due to the relative velocity between comoving 
observers on either side of the domain wall. This implies that for an observer looking across the domain wall the 
relation between the red-shift, z, and the scale factor, a, has to be modified to 

l + z = --2<P+. (16) 
a 

These effects are a distinguishing characteristic of these type of models, and could conceivably help distinguish 
them from quintessence-type models, for example. 



V. RESULTS AND DISCUSSION 



We have verified the accuracy of our formalism by computing the luminosity distance as a function of the red-shift 
of the source in two distinct ways. In the first approach (case I, dashed line in fig. 1) we assume that the background 
universe has cosmological parameters f2 A , QJ^ and H.q. On top of this we introduce a perturbation in the vacuum 
energy density parametrised by Af2° and calculate f(z) — TCocIl(z) /c using eqn. (|l3|). In the alternative ap proac h 
(case II, solid line in fig. 1) we assume that the cosmological parameters are fi^f , il^and H + AH (see eqns. (1(11)) 
with no perturbation. 

We have done the calculation for two distinct cosmological scenarios. Model A has fi^ = 1, 0° = and Af2° = 0.5 
while model B has f2^ = 0.3, 09 = 0.7 and Af2^ = 0.2. We can clearly see that the results obtained in either case 
are in very good agreement for both models. In fig. 2 we show the dependence of this agreement on the value of 
$0 = &{vo) for model B (here we took A^a > 0). We see that for |$q| "C 1 the relative agreement between the two 
methods for calculating the luminosity distance 



4(oo) 

is nearly proportional to |$o| 2 , as expected. 



(17) 
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Havingtested our method for calculating the luminosity distance we applied it in the context of the model described 
in Sect. |ll| and used type I supernovae in order to constrain the red-shift of the domain wall, z», and the fluctuation 
in the values of the cosmological constant in the outer region (parametrized by A£!° ). The high-redshift supernovae 
dataset of the Supernovae Cosmology Project was fit to the FRW magnitude redshift relation: 

mf = M B + 5\ogD L (z), (18) 

where -Dl(z) = H.qcIl(z) , Mb is the "Hubble-constant-free" £>-band absolute magnitude at maximum of a supernovae 
with a stretch factor s = 1, and m c g is the effective rest- frame B magnitude corrected for the width- luminosity relation. 
We assumed that fi m = 0.3 and fl\ = 0.7 in the inner region and we took Mb = —3.4 neglecting the uncertainty 
associated with the supernovae absolute magnitude calibration. (We have also checked that this particular choice 
does not affect our main results.) We have estimated the parameters z* and A£l^ using a x 2 statistical analysis. 

Fig. 3 displays our results in the z* versus Af2° plane. Note that for values of |<5| > 1 our approach based on 
linear perturbation theory ceases to be valid. The set of parameters for which this happens is denoted as region I. 
The region of parameter space that is allowed by the supernovae data (at 95% confidence) is denoted as region II, 
while the observationally excluded region is denoted as region III. 

We clearly see that, as expected, if z* is small only a small value of Afi^ is allowed. However, as the value of 
z* increases, the constraints on the values of AS!^ are much weaker, and there are essentially no constraints beyond 
z ~ 0.7. This happens essentially because the importance of the cosmological constant is smaller in the past than at 
the present time. 

We also see that negative values of Af2^ are excluded until significantly larger redshifts (z ~ 0.7) than positive ones 
(for which there are no constraints beyond z ~ 0.5). In fact the best fit model to the supernovae data has z* ~ 0.6 
and A£!° ~ 0.4. This is a significant result — it is an alternative (and perhaps intuitively clearer) way of saying that 
the high red-shift supernovae data does favour a relatively large positive cosmological constant. 

VI. CONCLUSIONS 

We have used type la supernovae data to constrain a simple toy model for spatial variations of the cosmological 
constant. Specifically, the model assumes that we are at the centre of a spherical region with a given set of cosmological 
parameters which is enclosed by a light domain wall and surrounded by another region where these parameters are 
different. This aims to be a very simple mimic model for cosmological models where the cosmological constant has 
space and/or time variations (such as, eg quintessence models). 

We have shown that the presently available supernovae data are only constraining at very low red-shifts. Never- 
theless, negative spatial variations (meaning a present-day value of the cosmological constant larger than the one at 
high red-shift) are much more constrained than positive variations. In fact, our best-fit model turns out to have a 
significantly larger vacuum energy density in the outer region than in the inner one. This is a clear indication that 
a positive vacuum energy density is favoured by the high red-shift data. Obviously larger and deeper datasets will 
significantly improve these constraints. 

Finally, we point out again that our analysis used a very simple toy model. While we do believe that the model still 
captures the crucial physics of the problem being discussed, it is clear that our analysis can be improved by resorting 
to a numerical simulation of the different domains. On the other hand, a careful discussion of the effects of these 
spatial variations of cosmological parameters in the CMB is also required. We shall return to these issues very shortly 
& 
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FIG. 1. Comparison between the luminosity distance, parametrized by f(z) = Hodt(z)/c, calculated in two distinct ways. 
In case I (dashed line) we assume a background universe with cosmological parameters £2°, Q,^ and Tto and a perturbation in 
the vacuum energy density parametrised by Afl^. In case II (solid line) we assume that the cosmological parameters are il* 
and H + AH with no perturbation. Model A has = 1, £2° = and AQ.% = 0.5 while model B has Q.° n = 0.3, 0° = 0.7 
and AD.% = 0.2. 
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FIG. 2. Dependence the relative agreement between the two methods for calculating the luminosity distances, on the value 
of | <E>o | for model B (here we took AS7a > 0). We see that e is nearly proportional to |$o| 2 , as expected. 
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